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ABSTRACT 


This  research  considers  communication  amongst  a 
community  of  earth  terminals  through  a  collection  of  low 
altitude  earth  oriented  satellites  (LASAT) .  In  one  proposed 
liASAT  system,  all  satellites  are  identical  and  simultaneous 
transmission  of  near  identical  signals  from  several 
satellites  creates  interference  at  the  earth  receivers. 
Frequency  hop  sequence  (FHS)  is  a  new  form  of  frequency 
hopping  (FH)  and  is  the  proposed  spread  spectrum  LASAT 
communications  technique  considered  in  this  research.  In 
military  applications,  the  perfomnance  of  the  communications 
system  in  the  presence  of  noise  only,  noise  and  jamming, 
noise  and  multipath  or  synchronization  errors  is  of 
interest.  FHS  offers  some  immunity  to  jamming  and 
multipath.  This  report  presents  the  calculated  noise 
performance  of  FHS  in  the  presence  of  jamming,  multipath  and 
synchronization  errors.  For  the  noise  only  case,  the 
calculated  probability  of  bit  error  of  a  FHS  system  using 
eight  symbols  for  a  10  dB  value  of  the  ratio  of  input  bit 
energy  to  noise  power  density  is  2  x  10'^. 
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INTRODUCTION 


Radio  communication  is  transitioning  from  analog  to 
digital.  The  quality  measure  of  a  digital  communication 
system  is  the  probability  of  error,  rather  than  output  SNR 
as  in  analog  communications.  Military  radio  communication 
is  concerned  with  the  effects  of  jamming.  This  has  produced 
the  spread  spectrum  in  which  transmitted  signals  are 
digital.  One  form  of  spread  spectrum  is  frequency  hopping 
(FH) . 

Techniques  and  performance  of  FH  are  well  documented 
[Ref.3].  FH  provides  immunity  to  interference  from  jamming 
and  multipath  effects  [Ref.4].  This  report  is  concerned 
with  the  application  of  a  particular  form  of  FH,  called 
Frequency  Hop  Sequence  (FHS) ,  to  a  particular  type  of 
military  communication  system  using  a  multitude  of 
low-altitude  satellites.  Chapter  II  gives  the  background 
about  FHS.  In  chapter  III,  the  performance  analysis  of  FHS 
is  provided.  Chapter  IV  contains  conclusions  and 
recommendations . 
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ZZ.  BACKGROUND 


In  a  LASAT  communication  system  proposed  by  Dr.  Glen  A. 
Myers,  satellites  provide  over  the  horizon  communication 
amongst  the  earth  stations  [Ref. 5].  This  research  considers 
that  proposal  in  which  all  satellites  are  identical  with  each 
consisting  of  a  single  wideband  transponder  (B  Hz  bandwidth) . 
Simultaneous  transmission  of  near  identical  signals  from  the 
several  satellites  creates  interference  at  the  earth  station 
receivers.  The  interference  is  similar  to  that  created  by 
multipath  signal  reception.  A  well  known  method  of 
eliminating  the  effects  of  multipath  signal  is  frequency 
hopping  (FH)  [Ref.4].  FH  also  provides  immunity  to  jamming  in 
a  hostile  environment  by  rec[uiring  the  jammer  to  spread  the 
transmitted  power  over  a  much  wider  band  (B  Hz) . 

A  new  form  of  FH,  called  Frequency  Hop  Sequence  (FHS) ,  is 
the  proposed  communication  technique  for  the  LASAT 
communication  system.  In  FHS,  a  group  of  bits  is  transmitted 
simultaneously  and  these  groups  or  symbols  are  distinguished 
during  transmission  by  assigning  a  unique  FH  sequence  to  each 
symbol.  In  all  forms  of  FH,  the  data  randomizes  the  hop 
pattern.  With  FHS,  the  data  mixes  the  sequences  rather  than 
simply  offseting  the  carrier  hop  by  hop.  The  notion  is  that 
it  may  be  more  difficult  to  recognize  and  hence  duplicate  an 
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entire  sequence  than  it  is  to  recognize  the  offset  of  a 
frequency  cell.  FHS  also  provides  immunity  to  interference  in 
the  form  of  multipath  and  jamming. 

The  problem  of  interest  in  this  report  is  the  noise 
performance  of  FHS  in  the  presence  of  interference  in  the  form 
of  multipath  and  jamming.  The  effect  of  synchronization  error 
is  also  considered. 
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SYSTEM  MODEL  AMD  RESULTS 


In  FHS,  a  group  of  k  bits  is  sent  simultaneously.  There 
are  then,  M=2''  possible  groups  or  symbols.  These  symbols 
are  distinguished  during  transmission  by  assigning  a  unique 
FH  sequence  to  each  symbol.  is  the  number  of  hops  per 
symbol,  is  the  number  of  frequency  cells  used,  b  is  the 
bandwidth  of  each  frequency  cell  and  B  is  the  total 
bandwidth  used  by  the  system.  All  of  these  system 
parameters  are  illustrated  in  Figure  1. 

The  receiver  of  the  system  has  2*^  parallel 
dehop/detector  subsystems  as  shown  in  Figure  2.  The  voltage 
of  interest  in  each  of  these  subsystems  is  represented  by  z„ 
where  m  =  1,2,3. ...M.  Each  of  these  subsystems  has  a 
synchonized  hopping  local  oscillator  (L.O.)  at  the  front  end 
with  a  unique  hopping  pattern.  If  the  hopping  pattern  of 
L.O. in  subsystem  one  is  the  same  as  the  hopping  pattern  of 
symbol  Sj,  then  the  output  of  the  intermediate  frequency 
(IF)  amplifier  in  subsystem  one  is  a  sinusoid  at  the  IF  with 
a  duration  T,  sec  and  amplitude  A  volts.  This  sinusoid  is 
envelope-detected  and,  to  improve  the  performance,  the 
voltage  at  the  output  of  the  envelope  detector  is  integrated 
for  a  symbol  duration  and  the  integrator  is  then  dumped. 
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Figure  1.  FHS  Frequency  versus  Time  Diagram 


When  the  resulting  output  voltages  of  M  subsystems  are 
compared  by  a  comparator,  the  decision  will  favor  subsystem 
one  or  Sj.  Then  a  symbol-to-bit  converter  produces  the 
desired  bit  stream. 

A.  NOISE  PERFORMANCE 

The  noise  performance  analysis  of  the  system  is  done  at 
the  output  of  the  envelope  detector.  It  is  assumed  that  the 
a  priori  probability  of  each  symbol  is  equal.  First,  the 
probability  of  making  a  symbol  error  is  found;  then,  with 
suitable  conversions,  the  probability  of  bit  error  is 
obtained. 

If  it  is  assumed  that  symbol  one  is  sent,  then  the 

probability  that  S,  is  detected  correctly  is  the  probability 

z,  is  greater  than  Zj  and  is  greater  than  Zj,  etc.  This 
can  be  expressed  as  Pgg=Pr  (  Zj  <  <  ZiD . Zm<ZiISi)  where 

n  is  the  and  symbol  and  Pr(A|B)  is  the  conditional 
probability  of  A  given  B.  Then  the  probability  of  making  an 

error  when  detecting  S,  is  Pk=1'-’ 

••*1  *1 

Pg^  =  l-JJ .  .  .  Jp(Zj^,Z2,2j. .  .ZM'Si)dZjdz2  —  dz,j  (1) 

0  0  0 

where  p(A,B|C)  is  the  joint  probability  of  A  and  B  given 
condition  C.  Since  all  subsystem  voltages  are  statistically 
independent,  the  conditional  probability  function  inside  the 
integrals  can  be  expressed  as, 
p(Zi,Zj, . ,z„!Si)  =  ft  P{z„!Si).  (2) 

IB-1 
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If  it  is  assimed  that  all  subsystems  except  the  one  with 


the  signal  will  have  the  same  statistical  behavior,  then 

Zi 

p{z„lSi)dZj,  are  identical  when  m  is  not  equal  to  one. 

0 

Thus  the  probability  of  making  an  error  when  detecting  Sj 
can  be  written  as, 

Pg,=J[l-[J’p(z„lSi)dzJ“-^]p(ZilSi)dZi  ,  m=2  or  3...or  M.  (3) 

0  0 

Now,  the  probability  density  function  of  the  voltage  at 
the  output  of  the  envelope  detector,  given  that  the  input  is 
zero  mean  Gaussian,  is  the  Rayleigh  density  function.  So, 
the  conditional  density  function  inside  the  parenthesis  can 
be  expressed  as, 

z  zi 

p(z_lS.)® — ^exp( - ,  m=2  or  3... or  M.  (4) 

2a 

where  Oi2®(Nq/2)  (2/Tg)  (2)  and  No/2  is  the  two-sided  power 
spectrum  of  the  noise  voltage  at  the  input  to  the  envelope 

detector.  The  integral  of  equation  (4)  is  evaluated  as 
*1 


Then,  the  expression  inside  the  parenthesis  in  equation  (3) 
can  be  written  using  the  binomial  expansion  as 

l-[l-exp(— exp(-m-^)  (6) 

20^2 


/p(2«ISi)dz„ 


/-^exp(— ;^)dz„=l-exp(— .  (5) 


where 


C 


M-l 

m 


(M-l)  I 
m!  (M-m-1) ! 
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Now,  p(Zi|Si)  is  a  Rician  density  function  with  the 
signal  amplitude  A.  So, 

"M-l 


^86=  /]C  ^  -1)  “*^exp  ( -m— )  -^exp-  (  ^  )  Iq  ( )  dz^  (7 ) 

J  ±r?  20j  «-  2a. ^  a.'^ 


Zi'  .  Z.  -  \T  ( 

OB-I  20^2  Oi 

Substituting  this  and  equation  (6)  into  equation  (3)  gives 

a2 


M-l 


exp[  -( 


Pae=E  (-1) 


in^l 


m 

m+1 


) 


] 


ni«l 


iri+l 


(8) 


Since  all  symbols  are  equally  likely,  equation  (8)  is  also 
the  probability  of  error  for  all  symbols. 

To  express  P,*  in  terms  of  the  ratio  of  input  symbol 
energy  to  noise  power  density,  define  the  symbol  energy  as 

Then,  the  ratio  of  input  symbol  energy  to  noise 


E,  = - - 


power  density  becomes  E3/No=A^/oi  .  When  this  expression  is 
substituted  in  equation  (8) ,  the  probability  of  symbol  error 

becomes 


M-l 


exp  { — 2- ) 

Ps8=e 

m-l 


(9) 


(m+1) 

since  a  symbol  consists  of  k  =  logj  M  bits,  then 
mistaking  one  symbol  for  another  can  cause  1  or  2  or  more  up 
to  k  bits  to  be  in  error.  In  fact,  the  number  of  bits  in 
error  is  just  the  Hamming  distance  of  the  two  symbols.  If 

all  symbol  errors  are  equally  likely,  - d-C is  the 

M-l  ^ 

average  number  of  bits  in  error  for  each  symbol  in  error. 

To  convert  the  probability  of  symbol  error  to  the 
probability  bit  error,  we  use  the  following  reasoning. 
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First,  the  transmission  of,  say,  1000  symbols  is  equivalent 
to  the  transmission  of  1000  x  k  bits.  Each  symbol  in  error 

-  k 

creates  an  average  number  - 52  bits  in  error.  So, 

M-1 

Using  the 


Pb,  becomes  x 

equation  (9)  for  gives 


•')  X  1 


M-l 


exp( — 

m+1 


m  V  EijlogjM 


2N„ 


(m+l) 


(10) 


where  we  have  let  E,=Efclog2  M.  Figure  3  is  plots  of  the 
probability  of  bit  error  versus  the  ratio  of  bit  energy  to 
noise  power  density  for  M  =  2,4,8,16.  For  the  two  symbol 
case,  the  receiver  was  simulated  and  the  probability  density 
functions  of  the  voltages  at  the  outputs  of  the  envelope 
detectors  for  signal  and  nonsignal  subsystems  were  obtained 
from  the  simulation  data.  These  probability  density 
functions  are  the  Rayleigh  density  function  for  the 
nonsignal  subsystem  and  the  Rician  density  function  for  the 
signal  subsystem  as  shown  in  Figure  4.  The  probability  of 
bit  error  for  a  particular  value  of  the  ratio  of  bit  energy 
to  noise  power  density  is  the  total  area  under  the  tails  of 
the  probability  density  functions  as  calculated  from 
equation  (10) . 

The  performance  curves  for  the  outputs  of  the  envelope 
detector  and  the  integrator  were  also  obtained  from  the 
simulation  data.  As  Figure  5  illustrates,  the  performance 
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Figure  3.  The  Probability  of  Bit  Error  versus  the  Ratio  of 
Bit  Energy  to  Noise  Power  Density  for 
M  =  2,4,8,16. 
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Figure  4.  The  Probability  Density  Functions  of  Voltage  at 
the  Output  of  the  Envelope  Detector  of  the 
Signal  and  Nonsignal  Subsystems.  M  =  2. 
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Figure  5 


Eb/NO 


The  Simulated  Values  of  Probability  of  Bit  Error 
versus  Bit  Energy  to  Noise  Power  Density  Ratio. 
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improvement  factor  of  the  integrator  is  about  3  dB  for  high 
ratios  of  bit  energy  to  noise  power  density,  but  it  is  less 
than  3  dB  for  low  ratios  of  bit  energy  to  noise  power 
density.  By  extrapolation,  this  same  performance 
improvement  can  be  obtained  for  the  4,8  and  16  symbol  cases. 
B.  THE  EFFECT  OF  SYNCHRONIZATION  ERROR 

If  the  hopping  L.O.  of  Figure  2  is  not  exactly  time 
synchronized  with  the  received  FHS,  then  the  output  of  the 

IF  amplifier  is  not  a  sinusoid  for  the  entire  symbol 
duration  (T,  seconds) .  So,  a  synchronization  error  of  e 
seconds  per  hop  (early  or  late)  causes  the  voltage 

amplitude  at  the  output  of  the  integrator  drop  from  AT, 

volts  to  BT,  volts  as  shown  in  Figure  6.  That  is,  the 

voltage  at  the  output  of  the  envelope  detector  drops  from  A 

volts  to  B=A(1 — — ^)  for  0<e<-^  .  Then,  the  energy  of  the 

symbol  at  the  output  of  the  integrator  decreases  from  E,  to 

A  €N 

E,=  (l — .  This  effectively  moves  the  probability  of 

bit  error  versus  the  ratio  of  bit  energy  to  noise  power 

density  curves  of  Figure  3  to  the  right  depending  on  the 

eNh 

value  of  the  ratio  e  =  — 2  . 

T. 

Figures  7,8,9,10  illustrate  the  probability  of  bit 
error  curves  for  a  possible  range  of  e  when  M  *  2,4,8,16 
respectively.  All  of  these  curves  should  be  moved  to  the 
left  in  accordance  with  Figure  5  to  include  the  improvement 
effect  of  the  integrator.  For  the  two  symbol  case,  the 
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Figure  6.  The  Voltage  Drop  at  the  Output  of 

the  Integrator  Due  to  Synchronization 
Error  in  the  Signal  Subsystem. 
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Figure  7.  The  Effect  of  Synchronization  Error  on 

System  Performance  for  M  =  2  and  for  Various 
Values  of  e. 
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Figure  8 


The  Effect  of  Sychronization  Error  on 
System  Performance  for  M  =  4  and  for  Various 
Values  of  e. 
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Figure  9.  The  Effect  of  Synchronization  Error  on 

System  Performance  for  :-i  =  8  and  for  Various 
Values  of  e. 


Figure  10.  The  Effect  of  Synchronization  Error  on 

System  Performance  for  M  =  16  and  for  Various 
Values  of  e. 
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receiver  with  a  synchronization  error  in  the  signal 
subsystem  was  simulated. 

The  probability  density  functions  of  voltage  at  the 
output  of  the  envelope  detector  for  the  signal  and  the 
nonsignal  subsystem  were  obtained  from  the  simulation  data 
and  are  shown  in  Figure  11.  The  Rician  density  function, 
which  is  the  density  function  in  the  signal  subsystem,  moves 
to  the  left  due  to  the  synchronization  error.  This 
increases  the  area  under  the  tails  and  so  also  the 
probability  of  error. 

C.  THE  EFFECT  OF  INTERFERENCE 

Interference  can  occur  because  of  multipath  effects  or 
jamming.  In  the  case  of  interference,  the  nonsignal 
subsystems  have  an  interference  signal  in  addition  to  the 
noise  input  to  the  system.  So,  the  probability  density 
functions  of  the  voltages  at  the  outputs  of  the  envelope 
detectors  become  the  Rician  density  functions  as  shown  in 
Figure  12  for  the  two  symbol  case.  The  probability  of  error 
is  the  area  under  the  tails  as  in  the  previous  cases.  Now, 
assume  that  the  effect  of  the  interference  is  to  increase 
the  mean  of  the  Rayleigh  distribution  at  the  output  of  the 
envelope  detector  of  the  nonsignal  subsystems.  The  Rician 
density  function  of  the  voltage  at  the  output  of  the 
envelope  detector  in  the  nonsignal  subsystem  is  then 
approximated  by  the  Rayleigh  density  function  such  that  the 
area  under  the  tails  of  the  probability  density  functions  or 
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Figure  11.  The  Probability  Density  Functions  of  Voltage  at 
the  Output  of  the  Envelope  Detector  of  the 
Signal  and  Nonsignal  Subsystems  in  Case  of 
Synchronization  Error.  M  =  2. 
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Figure  12.  The  Probability  Density  Functions  of  Voltage  at 
the  Output  of  the  Envelope  Detector  of  the 
Signal  and  Nonsignal  Subsystems  in  Case  of 
Interference.  M  =  2. 
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the  probability  of  error  does  not  change  as  shown  in  Figure 
13.  With  these  approximations,  the  analysis  can  be  done 
using  the  equations  developed  in  section  A  of  this  chapter. 
Two  cases  of  interference  are  considered:  interference  in 
all  nonsignal  subsystems  and  interference  in  only  one 
nonsignal  subsystem. 

1.  Interference  in  all  Nonsignal  Subsystems 

When  there  is  no  interference,  the  probability 
density  function  at  the  output  of  the  envelope  detector  in 
all  nonsignal  subsystems  is  the  Rayleigh  density  function 
with  a  mean  of  ^=Vii/20i  where  is  the  variance  of  the 

input  Gaussian  noise.  In  the  case  of  interference,  the  mean 

of  the  Rayleigh  distribution  at  the  output  of  the  envelope 
detector  of  the  interfered  subsystems  is  defined  as  y=z^ 

where  11  is  the  additional  dc  voltage  at  the  output  of  the 
envelope  detector  created  by  the  interference.  Since  the 

mean  of  the  Rayleigh  distribution  is  related  to  its  standard 
deviation  of  the  input  noise  by  y=^n/2o2  ,  then  the  same 
performance  analysis  approach  of  section  A  can  be  used  by 
taking  different  variances  for  the  signal  subsystem  (0^^) 

and  for  the  nonsignal  subsystems  (Oj^)  when  interference 
occurs.  The  equation  (4)  can  now  be  rewritten  as, 
z  2^ 

P(z„!Si)  =— ^exp( - ,  m=2,3,4, . M  (11) 

202^ 

where  y=Ui+<yi  and  since  y=^+ir  ,  ^=>/n/20i 

ir=v/n/20j  . 
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Figure  13.  The  Probability  Density  Functions  of  Voltage  at 
the  Output  of  the  Envelope  Detector  in  Case  of 
Interference  When  the  Rician  Density  Function 
is  approximated  by  the  Rayleigh  Density.  M  =  2. 
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(12) 


If  the  integral  of  the  equation  (11)  is  evaluated, 


”4 

Tp  ( z„  I  Si)  dz„=  T-^exp  ( )  dz^=l-exp  ( 


0  0  ®2  *''2 
then  the  equation  (7)  can  be  rewritten  as 


-z 


202^ 


■m-1 


p..=/E  ^  ^  ^  )  lo  ( -^ )  dzi .  (13 ) 

0“^  2o2‘‘  Oi  20i‘‘ 

When  this  integral  is  rearranged  and  evaluated,  the 


zf+A^ 


probability  of  symbol  error  becomes 


P»e-E  Cm'*'"  (-1) 

m*! 


I 


7ai= 


m-*-02^/oi‘ 


]  exp  t 


-  (m  A^/20i^) 
(m+Oj^/Oj^) 


1  . 


(14) 


Substituting  Eg/Nj,  =  A^/Oi^  ,  Eg  =  E^jlogjM  and  defining 

-  It 

z  =  and  including  the  factor  L=-; - nrE'^^n’^  * 

(M-1)  K  jfri 


probability  of  bit  error  becomes 


M-1 


Pb,*(L)  j;  (-1)“^^  [-5_]  expt  3  . 

n'+z  (m+z)2N(, 


(15) 


To  obtain  the  probability  of  bit  error  as  a 


function  of  the  ratio  of  bit  energy  to  noise  power  density 

and  signal  to  interference  ratio,  (SIR) ,  the  following 

a2  a 

relationships  are  defined:  SIR= — ,  a2=Oi+OT  ,  0^  =  ■  -  - —  • 

Using  these  relationships,  a  new  variable  r  is  defined  as 
(Eg/NQ)/SIR  .  Then  z  is  found  to  be  related  to  r  by  the 


following  quadratic  equation:  z=l+2>/r+r  . 

Figure  14  is  curve  of  z  versus  r  for  the  range 

of  interest  of  r.  Figures  15,16,17,18  show  the  family  of 

probability  of  bit  error  versus  the  ratio  of  bit  energy  to 

noise  power  density  curves  for  different  values  of  r. 
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Figure  14.  z=o\/q\  versus  r  =  (E./N^) /SIR 
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Figure  15.  The  Effect  of  Interference  on  System 

Performance  when  All  Nonsignal  Subsystems  have 
Interference  for  M  =  2.  (E,/Nj,)/SIR  is  Marked 
on  each  Curve. 
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Figure  16.  The  Effect  of  Interference  on  System 

Performance  when  All  Nonsignal  Subsystems  have 
Interference  for  M  =  4.  (E^/Nq) /SIR  is  Marked 

on  each  Curve. 
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Figure  17.  The  Effect  of  Interference  on  System 

Performance  when  All  Nonsignal  Subsystems  have 
Interference  for  M  *=  8.  (E^/Nq)/SIR  is  Marked 
on  each  Curve. 
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Figure  18.  The  Effect  of  Interference  on  System 

Performance  when  All  Nonsignal  Subsystems  have 
Interference  for  M  *  16.  (E,/Nq)/SIR  is  Marked 
on  each  Curve. 
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All  of  these  curves  should  be  moved  to  the  left  in 
accordance  with  Figure  5  to  include  the  performance 
improvement  effect  of  the  integrator. 

Once  z  is  determined  from  Figure  14,  the 
probability  of  bit  error  versus  the  ratio  of  bit  energy  to 
noise  power  density  curve  can  be  obtained  by  equation  (15) . 
In  section  D,  the  use  of  Figure  14  for  determining  the 
probability  of  bit  error  will  be  clarified  with  an  example. 
2.  Interference  in  One  Noneignel  Subsystem 
When  only  one  nonsignal  subsystem  has 
interference,  the  equation  (3)  of  section  A,  assuming 
subsystem  *'  '  had  interference,  can  be  modified  as 

^86=1-/  [/P(z„lSi)dzJ«-2  [|p(z2lSi)dz2]  P(ZilSi)dZi  (16) 

0  0  0 

where  m  =  3  or  4... or  M.  Then  equation  (8)  of  section  A  can 
be  rewritten  as, 

P8e=l-t/-^exp-(^i^)I(,(^dzJ  (fi(m,Zi) -f2(m,Zi) )  (17) 

<j  Cl  2ci  Oj 

where 

M-2 

f  1  (m,  Zi)  =52  (-l)®exp-( - i)  and 

B-o  20i 

f2(m,z,)=V  (-i)«  exp[-(-5L+-%)-^]  . 

B-O  Oi  o|  2 

The  probability  of  symbol  error,  after 

evaluating  the  integral  and  carrying  out  the  multiplication 

in  equation  (17),  has  the  foirm 

Pg.*!- (fi  (m) -fj  (m) )  (18) 
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where 


exp[  -( 


m 


gi(ni)=52  C/-2  (-1)“ 

IB'O 

M-2 

g^Cm)  =5;c“-"(-l)"  [■ 


a2 


m+1  20^2 


m+l 


and 


0  exp[-^(. 


)] 


(m+l) +Oi/o|  20i  (m+l)+of/o| 

After  carrying  out  the  substitutions  for  the 

symbol  energy  the  probability  of  bit  error  as  a  function  of 

the  ratio  of  bit  energy  to  noise  power  density  and 
2/2 

z^Oa/oj  IS 

k 


Pb.-(  (”>))! 


(19) 


where 


gi(m)=53  (-1)“ 


M-2 

=E 

IB>0 

M-2 

=E 

IB«0 


exp[ 

m+l  2Nq 
m+l 


and 


g2(m)=^cr"  (-!)•“  exp  f-7- )  1  . 

^  “  (m+l) +1/2  2No  (m+l) +1/2 


Figures  19,20,21,22  show  the  family  of 
probability  of  bit  error  versus  the  ratio  of  bit  energy  to 
noise  power  density  for  different  values  of  r  =  E,/No/SIR 
when  M  =  2,4,8,16  respectively.  All  of  these  curves  should 
be  moved  to  the  left  in  accordance  with  Figure  5  to  include 
the  improvement  effect  of  the  integrator.  As  before, 
equation  (19)  should  be  used  by  determining  z  from  Figure 
14. 
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Figure  19.  The  Effect  of  Interference  on  System 

Performance  when  One  Nonsignal  Subsystem  has 
Interference  for  M  =  2.  r*(E,/No) /SIR  is  Marked 


on  each  Curve. 


Figure  20.  The  Effect  of  Interference  on  System 

Performance  when  One  Nonsignal  Subsystem  has 
Interference  for  M  =  4,  r=(E,/No) /SIR  is  Marked 

on  each  Curve. 


Figure  21.  The  Effect  of  Interference  on  System 

Performance  when  One  Nonsignal  Subsystem  has 
Interference  for  M  =  8.  r=  (E,/No) /SIR  is  Marked 
on  each  Curve. 
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Figure  22.  The  Effect  of  Interference  on  System 


Performance  when  One  Nonsignal  Subsystem  has 
Interference  for  M  =  16.  r=(E,/No) /SIR  is  Marked 

on  each  Curve. 
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D.  EXAMPLES 


The  use  of  equations  for  the  synchronization  error  and 
the  interference  cases  are  illustrated  with  the  following 
examples. 

The  system  parameters  used  in  the  examples  are: 
bit  rate  =  2000  bits/sec. 
number  M  of  symbols  =  4 
number  (N,,)  of  hops  per  symbol  =  10 

Example  1:  Determine  the  probability  of  bit  error  for  a 
10  dB  ratio  of  bit  energy  to  noise  power  density  in  the  case 
of  25  micro  seconds  synchronization  error  and  no 
interference. 

To  find  the  probability  of  bit  error,  the  parameter  e  needs 
to  be  calculated.  From  the  system  parameters,  e  is  0.25. 
Then  Eb/No®  (1-e)  ^  Eb/N(,  is  substituted  in  equation  (10)  and 
the  probability  of  bit  error  is  0.06.  Including  the  effect 

of  the  integrator  in  accordance  with  Figure  5,  the 
probability  of  bit  error  becomes  0.01. 

Example  2:  Determine  the  probability  of  bit  error  for 
a  14  dB  ratio  of  bit  energy  to  noise  power  density  if  all 
nonsignal  subsystems  had  interference  with  r  equal  -5  dB 
and  there  was  no  synchronization  error. 

To  find  the  probability  of  bit  error,  the  parameter  z  as 
read  from  Figure  14,  which  is  1.7.  Then,  substituting  this 
value  of  z,  the  value  of  the  ratio  of  bit  energy  to  noise 
power  density,  the  values  of  k  and  M  in  equation  (15) ,  the 
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the  integrator  in  accordance  with  Figure  5,  the  probability 
of  bit  error  becomes  5  x  10'*. 

Instead  of  directly  using  the  equations,  the  same 
results  could  also  be  obtained  by  using  the  probability  of 
bit  error  curves  in  both  cases  since  the  curves  for  e  =  0.25 
and  r  =  -10  dB  are  provided  as  Figures  8  and  16 
respectively. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

System  performance  improves  as  the  number  of  subsystems  or 
symbols  increases,  at  the  expense  of  more  hardware.  The 
performance  improvement  due  to  an  increase  in  the  number  of 
symbols  is  larger  in  value  for  larger  values  of  the  ratio  of 
bit  energy  to  noise  power  density.  The  integrator  in  the 
receiver  improves  the  performance  by  about  3  dB  for  large 
values  of  the  ratio  of  bit  energy  to  noise  power  density,  but 
the  improvement  gets  less  than  3  dB  for  small  values. 

The  synchronization  error  degrades  the  system  performance 
regardless  of  the  nvimber  of  symbols  the  system  uses.  For  the 
same  amount  of  synchronization  error,  the  performance 
degradation  worsens  as  the  niimber  of  hops  per  symbol 
increases,  as  expected.  If  large  synchronization  errors  are 
expected,  the  number  of  hops  per  symbol  should  be  chosen 
carefully. 

Interference  in  all  nonsignal  subsystems  or  in  one 
subsystem,  degrades  system  performance.  When  only  one 
nonsignal  subsystem  has  interference,  the  degradation  worsens 
as  the  number  of  subsystems  decreases.  For  the  same  amount  of 
interference  power,  the  interference  in  all  nonsignal 
subsystems  degrades  the  system  performance  more  than  the 
interference  in  only  one  nonsignal  subsystem. 


39 


The  composite  effect  of  interference  and  synchronization 
error  is  to  further  degrade  the  system  performance.  When  one 
nonsignal  subsystem  has  interference  and  there  is 
synchronization  error,  both  the  number  of  hops  per  symbol  and 
the  number  of  subsystems  that  the  system  uses  become  important 
design  parameters. 

B.  RECOMMENDATIONS 

It  is  recommended  that  the  results  obtained  for  the 
performance  of  the  FHS  system  considered  in  this  research  be 
verified  by  a  hardware  realization. 
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APPENDIX 


MATLAB  SIMULATION  CODE 

%This  program  simulates  the  receiver  of  the  FHS  system  for  % 
%two  symbol  case.  For  a  given  input  bit-energy-to-noise  % 
%density  ratio,  the  output  of  the  program  is  the  % 

%probability  of  bit  error  for  that  bit-energy-to-noise  % 

%density  ratio.  % 

function  nO  =  noisevec(c) 

%Forms  the  band  limited  white  noise  vector  with  variance  10% 
rand ('normal')  ; 

M«rand(c,4000)  ; 

B=(sqrt (10) /std(A) )  *  M  ; 
f=[0  .01  .02  .03  .04  .05.. 

.06  .07  .08  .09  .1  .2 

.3  .4  .5  .6  .7  .8  .9  1]  ; 

k=[l  1111111111.. 

000000000] 

[b,a]=yulewalk(15,f ,k)  ; 

C=f ilter (b,a,B)  ; 

n0=reshape( (sqrt (10) /std(C) )  *  C, 20, 200); 

function  F=sigenv(m)  ; 

A=sqrt(10*m)  ;%Form  the  Eb/NO 
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for  n=l:20 


t 


nO=noisevec(l)  ;%Get  the  noise  vector 

%ForTn  the  envelope  of  signal+noise% 
Es=sqrt((A+nO(:, 1:100)) . *2+(A+nO ( : , 101: 200) ) .^2) ] 
%Integrate  for  a  symbol  duration  and  then  dump% 
Is(n, : )=sum(Es) . /20  ; 

end  ; 

F=reshape(Is, 2000,1)  ; 

function  E=noiseenv(N)  ; 

for  d=l:20  ; 

n0=noisevec(l)  ; 

%Form  the  envelope  of  the  noise% 

En=sqrt(n0(: ,1:100) .^2+n0(: ,101:200) .^2) ; 
%Integrate  for  a  symbol  duration  and  dump% 

In(d, :)=sum(En) ./20 
end 

E=reshape (In, 2000, 1) 

%  Main  Program 
F=sigenv(m) 

E=noiseenv (20) 
ms=mean(F) 
mn=mean(E) 

Vt= (ms+mn) /2 
i=length ( f ind ( F<Vt ) ) 
j=length(f ind(E>Vt) ) 

Pm=i/2000 


%  Signal+noise  vector 
%  Noise  only  vector 
%  Find  the  mean  of  signal+noise 
%  Find  the  mean  of  noise  only 
%  Find  the  ML  threshold 

%  Find  probability  of  miss 
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probability  of  false  alarm 
probability  of  error 
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